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The optimal precision of frequency measurements in the 
presence of decoherence is discussed. We analyse different 
preparations of n two-level systems as well as different mea- 
surements procedures. We show that standard Ramsey spec- 
troscopy on uncorrelated atoms and optimal measurements on 
maximally entangled states provide the same resolution. The 
best resolution is achieved using partially entangled prepara- 
tions with a high degree of symmetry. 



The rapid development of laser cooling and trapping 
techniques has opened up new perspectives in high pre- 
cision spectroscopy. Frequency standards based on laser 
cooled ions are expected to achieve accuracies of the or- 
der of 1 part in 10 14 — 10 18 Q. In this letter we discuss 
the limits to the maximum precision achievable in the 
spectroscopy of n two level atoms in the presence of de- 
coherence. This question is particularly timely in view 
of current efforts to improve high precision spectroscopy 
by means of quantum entanglement. 
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FIG. 1. Schematic representation of Ramsey spectroscopy 
with uncorrelated particles. 

In the present context standard Ramsey spectroscopy 
refers to the situation schematically depicted in Fig. 1. 
An ion trap is loaded with n ions initially prepared in 
the same internal state |0). A Ramsey pulse of frequency 



u> is applied to all ions. The pulse shape and duration 
are carefully chosen so that it drives the atomic transition 
|0) <-> |1) of natural frequency u>q and prepares an equally 
weighted superposition of the two internal states |0) and 
|1) for each ion. Next the system evolves freely for a 
time t followed by the second Ramsey pulse. Finally, the 
internal state of each particle is measured. Provided that 
the duration of the Ramsey pulses is much smaller than 
the free evolution time t, the probability that an ion is 
found in |1) is given by 



P = (l + cosA*)/2. 



(1) 



Here A = u> — u>q denotes the detuning between the clas- 
sical driving field and the atomic transition. 

This basic scheme is repeated yielding a total duration 
T of the experiment. The aim is to estimate A as ac- 
curately as possible for a given T and a given number 
of ions n. The two quantities T and n are the physical 
resources we consider when comparing the performance 
of different schemes. The statistical fluctuations associ- 
ated with a finite sample yield an uncertainty AP in the 
estimated value of P given by 



AP = y'Pil- P)/N 



(2) 



where N = nT/t denotes the actual number of exper- 
imental data (we assume that N is large). Hence the 
uncertainty in the estimated value of ujq is given by: 
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(3) 



This value is often referred to as the shot noise limit [0 . 

The theoretical possibility of overcoming this limit has 
been put forward recently Jj^EJ]. The basic idea is to 
prepare the ions initially in an entangled state, which for 
small n seems to be practical in the near future. To see 
the advantage of this approach, let us consider the case 
of two ions prepared in the maximally entangled state la] 
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This state can be generated, for example, by the initial 
part of the network illustrated in Fig. 2. A Ramsey pulse 



1 



on the first ion is followed by a controlled-NOT gate || . 
After a free evolution period of time t the state of the 
composite system in the interaction picture rotating at 
the driving frequency ui reads 



(100) 



-2iAt 



|H»A/2. 



(5) 



The second part of the network allows to disentangle the 
ions after the free evolution period. The population in 
state |1) of the first ion will now oscillate at a frequency 
2A 



P 2 = (1 + cos2Ai)/2. 



(6) 
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FIG. 2. Spectroscopy with two maximally entangled parti- 
cles. The particles are entangled and disentangled by means 
of "controlled-NOT" gates |j. 

This scheme can be easily generalised to the n ion case 
by a sequence of controlled-NOT gates linking the first 
ion with each of the remaining ones. In this way, a max- 
imally entangled preparation of n ions of the form 



I*) 
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(7) 



is generated. The final measurement on the first ion, after 
the free evolution period and the second set of controlled- 
NOT gates, gives the signal 



P n = (1 + cosnAt)/2 



(8) 



The advantage of this scheme is that the oscillation fre- 
quency of the signal is now amplified by a factor n with 
respect to the case of uncorrelated ions and the corre- 
sponding frequency uncertainty is 



\Suj \ 



1 



n-JTt 



(9) 



Note that this result represents an improvement of a fac- 
tor 1 / \Jn over the shot noise limit (0) by using the same 
number of ions n and the same total duration of the ex- 
periment T 0] and it was argued that this is the best 
precision possible ||. 

Let us now examine the same situation in a realistic 
experimental scenario, where decoherence effects are in- 
evitably present. The main type of decoherence in an 
ion trap is dephasing due to processes that cause ran- 
dom phase changes while preserving the population in 
the atomic levels. Important mechanisms that result in 
dephasing effects are collisions, stray fields and laser in- 
stabilities. We model the time evolution of the reduced 



density operator for a single ion p in the presence of de- 
coherence by the following master equation f|: 

p(t) = iA (p\l)(l\ - \l)(l\p) + 7 (a z pa z - p) (10) 

the density matrix decay [Actually, our analysis is not 
restricted to this particular model but holds for any 
process where off-diagonal elements decay exponentially 
with time.] Equation ( [To| ) is written in a frame rotating 
at the frequency u>. By <j z = |0)(0| — |1)(1| we denote a 
Pauli spin operators. Here we have introduced the decay 
rate 7 = l/rd ec , where Td ec is the decoherence time. For 
the case of standard Ramsey spectroscopy this will give 
rise to a broadening of signal (1): 



P = (1 + cosAie" 7 *)/2. 



(11) 



As a consequence the corresponding uncertainty in the 
atomic frequency is no longer A-independent. We now 
have 
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(12) 



In order to obtain the best precision it is necessary to 
optimise this expression as a function of the duration of 
each single measurement t. The minimal value is attained 
for 



At = kir/2 (k odd) 

t = Tdec/2 



(13) 



provided that T > T,j ec /2. Thus the minimum frequency 
uncertainty reads 



1^0 1 opt = V — 
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(14) 
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FIG. 3. Frequency uncertainty \8luq\ as a function of the 
duration of a single shot t for maximally entangled and un- 
correlated particles. Note that the minimum uncertainty is 
exactly the same for both configurations. 
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For maximally entangled preparation the signal (||) in 
the presence of dephasing is modified as follows: 

P n = (1 + cos jiA ie-™ 7 *)/2 (15) 

and the resulting uncertainty for the estimated value of 
the atomic frequency is now minimal when 



At = kn/2n {k odd) 
t = T dec /2n 



(16) 



Interestingly, we recover exactly the same minimal uncer- 
tainty as for standard Ramsey spectroscopy (14). This 
effect is illustrated in Fig. 3. The modulus of the fre- 
quency uncertainty \5u>o\ is plotted as a function of the 
duration of each single experiment t for standard Ramsey 
spectroscopy with n uncorrelated particles and for a max- 
imally entangled state with n particles. In the presence of 
decoherence both preparations reach the same precision. 
This result can be intuitively understood by considering 
that maximally entangled states are much more fragile 
in the presence of decoherence: their decoherence time is 
reduced by a factor n and therefore the duration of each 
single experiment t has also to be reduced by the same 
amount. Moreover, the limit (^) represents the best ac- 
curacy for both uncorrelated and maximally entangled 
preparations and cannot be overcome by engineering dif- 
ferent kinds of measurements as will be shown below. 

Note that the problem addressed in precision spec- 
troscopy (i.e. the measurement of small atomic phase 
shifts) maps onto that of statistical distinguishability of 
nearby states, analyzed by Wootters (l(J and generalized 
by Braunstein and Caves ]ll| , |i2| ]. They have provided 
an upper bound for the precision in the estimation of 
a given variable that parametrizes a family of quantum 
states. In our case this variable is the detuning A. More- 
over, the optimal measurements always correspond to a 
set of orthogonal projectors in the n ions Hilbert space. 
It is worthwhile pointing out the generality of this re- 
sult in the sense that it accounts for any possible joint 
measurement on the n particles and any method of data 
analysis. When the Braunstein and Caves optimization 
procedure is applied to either uncorrelated or maximally 
entangled preparations of n ions it yields the same limit 

However, we will show in the following that with cer- 
tain partially entangled preparations one can overcome 
the limit (|l4|). Let us analyze first the case of general- 
ized Ramsey spectroscopy, namely a scheme where the 
operator 



fc=i 



(17) 



is measured after the free evolution period. In ( JT7| ) a x = 
|1)(0| + |0)(1| denotes the Pauli spin operator and the 
superscript k refers to the fcth particle. 

We can easily evaluate the expectation values of the 
operators S x and S x in terms of the corresponding quan- 
tities in the absence of decoherence: 



(18) 



(S 2 x )=n+( 
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)=n + e-^{{Sl{ 1 = 0))-n). 



(19) 



Finally, the resulting uncertainty in the atomic fre- 
quency is given by 
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where (AS X ) 2 = (S 2 ) - (S x ) 2 and N denotes the total 
number of measurements performed during the total time 
T . A straightforward calculation leads to an optimal du- 
ration of each measurement t opt given by the solution of 
the following equation [O 



n[l + (2 7 t - l)e 2 ^] = (AS v (t = 0)) 2 
and the corresponding sensitivity takes the form 



\Slo 



opt 



2n-ye 2 "< t ° 



T(S x (t = 0)) ; 



(21) 



(22) 



where (S x (t = 0)) and (S y (t = 0)) 2 refer to the initial 
state preparation. The subscript opt emphasizes that 
optimization with respect to At has already been taken 
into account, the minimum value being achieved for At = 
Tt/2. Notice that each initial preparation is optimized 
for different values of the single shot time. We can now 
state a lower bound for the precision attainable within 
this approach as follows: 



\8u> 



Opt 



> 



2 7i7 



T(S x (0))' 



> 



nr dec T 



(23) 



Compared with the results above, a maximum improve- 
ment of 1/y/e in the resolution can be achieved. Thus 
we have found that the bound (Il4)can be overcome for 
certain partially entangled states~[ fjj . 

We now analyze the best precision that can be achieved 
when optimizing the experiment with respect to both the 
initial state preparation and the final measurement. The 
problem has been studied by means of a numerical op- 
timization procedure and we have restricted ourselves to 
small numbers of ions in the trap. The initial state prepa- 
ration which leads to the best precision is of the form 
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(24) 



k=0 



where \k) denotes an equally weighted superposition of 
all states of n ions which contain cither a number A; or a 
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number n — k of excited states. By [_• • -J we denote the 
corresponding integer part. The coefficients au can be 
chosen to be real. For example, the corresponding family 
of states for n = 4 reads 

\ip4) = a (|0000> + |1111)) + oi(|0001) + |0010) (25) 
+ |0100) + 1 1000) + |1110) + 1 1101) + |1011) + |0111)) 
+a 2 (|0011) + |0101) + 1 1001) + |1100) + |1010) + |0110» . 
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FIG. 4. The optimum percentual improvement in the pre- 
cision relative to the limit ( |l^ ) as a function of the number 
of ions n. Solid line: Numerical optimization with respect 
to the initial preparation and application of the Braunstein 
and Caves algorithm for determinig the optimal measurement. 
Dashed line: Optimized initial preparation and generalized 
Ramsey spectroscopy as the final measurement. 

Note that this family of states exhibits a high degree 
of symmetry: it is completely symmetric under permuta- 
tions of the n ions and under exchange of the excited and 
the ground state for each ion. The optimum percentual 
improvement in the precision relative to the limit jli| ) as 
a function of the number of ions n is shown in Fig. 4. 
The solid curve shows the improvement obtained by op- 
timizing both the inital preparation and the final mea- 
surement using the algorithm of Braunstein and Caves. 
The dashed line exhibits the improvement obtained by 
optimizing only the initial preparation and performing 
the measurement given in Eq. ( |l7| ) corresponding to gen- 
eralized Ramsey spectroscopy. The improvement ob- 
tained by optimising the measurement is rather small. 
The question whether the (dashed) curve corresponding 
to Ramsey spectroscopy asymptotically saturates at the 
theoretical limit (^J) or below remains to be addressed. 
Moreover, whether the curve corresponding to Braun- 
stein and Caves optimization saturates at the same value 
or higher than the Ramsey curve is an open question. 

In conclusion, we can state that standard Ramsey spec- 
troscopy is optimal for uncorrclated particles both in 
the presence and in the absence of decoherence effects. 
The use of maximally entangled states does not provide 



higher resolution as compared to using independent par- 
ticles when decoherence is present. The best sensitivity 
is achieved when the ions are initially prepared in highly 
symmetric but only partially entangled states. 
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